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Abstract 

In  this  study  we  apply  an  interface  sensitive  ellipsometry 
technique  to  study  the  evolution  of  the  Si-SiOj  interface  as  a 
function  of  high  temperature  annealing  (750‘’-1100“C)  .  Essentially, 
the  ellipsometry  technique  embodies  the  use  of  liquids  that 
refractive  index  match  with  the  bulk  film  thereby  removing  the 
optical  response  of  the  overlayer  and  greatly  enhancing  sensitivity 
to  the  interface.  According  to  both  time  and  temperature  of 
anneal,  distinct  modes  of  behavior  are  observed  for  the  evolution 
of  the  interface.  For  short  anneal  times  a  rapid  change  in  the 
interface  is  observed  that  correlates  with  the  disappearance  of 
protrusions,  followed  by  a  slower  change  that  correlates  with  the 
disappearance  of  the  suboxide.  At  high  temperatures  viscous 
relaxation  dominates  while  at  low  temperatures  the  suboxide 
reduction  is  apparent.  A  model  for  the  interface  in  terms  of 
chemical  and  physical  interface  processes  is  proposed  and  model 
parameters  are  compared  with  literature  results. 
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Introduction 


The  Si/SiOj  interface  continues  to  be  a  topic  of  intensive 
study.  Of  particular  importance  are  properties  of  the  interface 
such  as  structure,  thickness,  uniformity  and  phase  composition 
which  have  been  explored  by  numerous  methods  including  transmission 
electron  microscopy,  TEM  (1-2),  infrared  spectroscopy  (3),  low- 
energy  electron  diffraction  (4)  ,  scanning  tunneling  microscopy  (5)  , 
ellipsometry  (6-10),  and  a  variety  of  surface  spectroscopies  (11). 
Presently,  ultra-thin  SiOj  films  (thickness  less  then  30  nm)  find 
extensive  application  in  submicron  integrated  circuit  technology 
where  even  a  small  degree  of  interfacial  microroughness  or 
nonuniformity  can  alter  device  performance  and  reliability.  In 
order  to  improve  interfacial  quality,  such  processing  procedures  as 
two-step  oxidation  process  (12-13)  and  post-oxidation  annealing 
(14-15)  have  been  proposed.  Both  methods  include  annealing  of  the 
grown  oxides  in  a  nonoxidizing  atmosphere  such  as  argon  or 
nitrogen.  High-resolution  TEM  was  used  to  investigate  the  influence 
of  the  annealing  on  interface  smoothness  (12-13,  16),  but  the 
results  are  ambiguous  and  on  some  points  contradictory.  Yet  there 
seems  to  be  general  agreement  that  annealing  improves  both 
electronic  properties  and  interface  flatness. 

Spectroscopic  ellipsometry  has  been  applied  to  study  interfaces 
under  transparent  films(17),  and  recently  we  reported  an  interface 
enhanced  immersion  spectroscopic  technique (18) ,  which  is  sensitive 
to  such  interface  characteristics  as  microroughness,  thickness  and 
phase  composition.  In  the  present  study  we  apply  this  new  method 
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to  determine  the  properties  and  reactions  at  the  Si/SiO,  interface 
during  the  post-oxidation  annealing. 

Experimental  Procedures  and  Data  Analysis 

Single-crystal  (100)  oriented  2  Dcm  p-type  silicon  wafers  were 
cleaned  using  a  slightly  modified  RCA  procedure  (19)  and  thermally 
oxidized  in  a  fused  silica  tube  furnace  in  clean  dry  oxygen  to 
about  24-27  nm  SiOj  at  SOO^C.  Upon  completion  of  a  particular 
oxidation  one  sample  was  removed  without  annealing  as  a  control  and 
the  others  were  annealed  in  a  clean  nitrogen  atmosphere.  We  used 
annealing  temperatures  and  times  in  the  range  750-1100“C  and  1-120 
min,  respectively. 

In  order  to  investigate  the  Si/Si02  interface  we  have  applied 
a  novel  enhanced  sensitivity  spectral  and  variable  angle  of 
incidence  immersion  ellipsometry  technique  (18) .  The  key  feature 
of  this  technique  is  that  spectroscopic  and  variable  angle  of 
incidence  ellipsometry  measurements  are  performed  in  a  transparent 
liquid  ambient  that  has  optical  properties  very  close  to  the  those 
of  the  SiOj  overlayer.  Therefore,  the  overlayer  film  is  "optically” 
(not  physically)  eliminated  and  the  probing  light  beam  becomes 
highly  sensitive  to  the  interface  properties. 

Generally  it  is  difficult  to  achieve  a  perfect  refractive 
index  match  for  the  liqur’d  ambient,  n„,  and  the  SiOj  overlayer,  n^,, 
over  a  broad  spectral  range.  Therefore  deviations  are  accounted 
for  in  the  analysis (18) . 

Carbon  tetrachloride  (CCI4)  is  a  suitable  immersion  liquid  for 


matching  index  to  SiOj  films.  The  refractive  index  of  the  liquid 
ambient,  n^,  was  calculated  from  a  Cauchy  dispersion  formula,  taking 
into  account  temperature: 


(1) 


Where  X  in  A  is  the  wavelength  of  probing  beam,  parameters  n*  = 
1.4427,  a  =  5.15x10^A-2  at  T=24.8“C  (20)  and  an/aT=0 . 00055  at  T=20°C 
for  pure  CCI4  (21)  . 

It  is  known  that  the  refractive  index  of  thin-film  SiOj 
over  layers,  n^,,  depends  on  thickness,  1,^,  oxidation  temperature, 
T„,  annealing  temperature,  T^,  and  time,  t„,  index  of  ^he  substrate 
orientation,  N,  and  preparation  conditions  such  as  preoxidation 
cleaning  and  oxidation  ambient.  The  spectral  dependence  of  the 
non-annealed  thin-film  siOj  refractive  index,  n^®(Lo»,  To,,X)  ,  for  a 
substrate  with  particular  orientation,  e.g  (100),  was  calculated 
from  a  single  term  Sellmeier  approximation: 


(2) 


where  A(Lo,,  T<„,  N)  and  \,(L^,  T^,,  N)  are  dispersion  parameters  that 
are  also  dependant  on  the  overlayer  thickness  and  preparation 
conditions  and  with  values  A=«1.15  and  X<,=92.3  nm  for  25-35  nm  thick 
thermally  grown  at  800®C  on  Si(lOO)  SiOj  films  (10). 

It  has  been  observed  that  the  refractive  index  of  the  SiOj 
over layer  undergoes  a  relaxation  during  annealing  (22-25) . 
Ellipsometric  measurements  of  the  thermal  relaxation  of  the  SiOj 
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single  film  effective  refractive  index,  i.e.  without  taking  the 
interface  into  consideration,  shows  (23)  that  two  relaxation 
processes  are  involved:  rapid  (-10-20  min.)  initial  relaxation  and 
slow  exponential  decay  to  the  refractive  index  of  fully  relaxed 
oxide.  The  slow  exponential  decay  is  because  of  the  elastic  strain 
relaxation  and  well  described  as: 


n 


OV 


t^,N,X)  =  n^d)  +  [riot  (A.)  -n^(l)]exp 


(3) 


where  n,  is  the  refractive  index  for  fully  relaxed  oxide,  and  the 
empirical  parameter  7j=0.17  (22),  which  is  a  measure  of  the  strength 
of  the  coupling  between  the  index  and  relaxation.  The  relaxation 
time  for  the  oxide  density  is  expressed  as: 


toSXp 


kT. 


an  / 


(4) 


where  T(,-2x10‘“  min  and  E,-6  eV  (22)  .  We  have  assumed  that  the 
refractive  index  of  the  fully  relaxed  oxide  film  is  equal  to  that 
for  the  bulk  oxide,  and  thus  the  spectral  dependence  nf(X)  was 
calculated  with  eqn. (2)  using  A=1.099  and  Xo=92.27  (10)  which 
corresponds  to  bulk  SiOj.  The  observed  initial  rapid  relaxation  of 
the  effective  refractive  index  (23)  can  be  attributed  to  the 
evolution  of  the  interface,  and  a  detailed  model  for  this  will  be 
discussed  below. 

In  order  to  obtain  unknown  interface  parameters,  we  used  a 
Marquardt  non-linear  best  fit  algorithm  which  minimizes  the  value 
of  the  error  function: 
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Q 


(5) 


*  •  J 


where  P  is  a  vector  of  unknown  interface  parameters,  is  the 
photon  energy,  0i  is  the  angle  of  incidence,  and  the  subscripts  cal 
and  exp  refer  to  calculated  and  experimentally  derived  values.  A"* 
and  are  the  values  obtained  using  the  vector  P  from  expanded 
Fresnel  formulas  (26)  and  a  matrix  algorithm  for  the  multilayer 
system  with  optical  parameters  for  the  CCI4  ambient  and  Si02 
overlayer  calculated  from  the  eqns. (1-4)  and  the  previously 
determined  thickness  of  the  SiOj  overlayer.  The  Bruggeman  effective 
medium  approximation  (27) ,  BEMA,  was  used  to  calculate  the 
effective  dielectric  function  of  the  interface,  which  was  modeled 
as  a  mixture  of  constituents  (discussed  later)  with  known  optical 
properties.  Volume  fractions  of  the  constituents  were  calculated 
from  assximptions  about  the  geometry  of  the  interface  that  are 
discussed  below,  and  the  geometrical  parameters  have  been  included 
as  unknown  parameters. 

A  commercially  available  vertical  ellipsometer  bench  was 
modified  to  become  a  rotating  analyzer  spectroscopic  ellipsometer 
(the  essential  features  were  previously  described  (28))  and 
calibrated  according  to  a  published  procedure  (29) .  A  specially 
designed  variable  angle  of  incidence  immersion  cell  that  fits  on 
the  ellipsometer  stage  was  used  (18) .  For  each  sample 
spectroscopic  ellipsometry  was  performed  both  jn  air  and  CCI4  in  the 
2. 4-4.0  eV  range  of  photon  energies,  E,  and  at  two  angles  of 
incidence,  <p,  of  12°  and  75°.  The  measurements  yield  the 
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ellipsometric  angles  A“‘’(E,0,T„,t„)  (Fig.l)  and  t„)  as  a 

function  of  annealing  conditions  (T„,  t„)  .  As  was  previously 

discussed  (18)  the  sensitivity  of  ^  to  interface  changes  is 
considerably  less  than  for  A,  hence  only  A  is  used  to  correlatee 
with  process  changes.  From  the  measurements  done  in  air,  the 
thickness  of  the  SiOj  over layer  was  obtained  without  significant 
error,  because  there  is  very  low  interface  sensitivity  for  the 
measxirements  in  air  (18)  . 

In  order  to  obtain  interface  dependencies,  which  characterize 
the  evolution  of  the  interface  via  annealing,  we  have  calculated  an 
effective  relative  interface  parameter: 

e„)  .  c„)  -  Af'  -  SA“^(r„,  c„)  («) 

where  A“’’(T^,  t„)  is  the  experimental  ellipsometric  angle  A  at  an 
annealing  temperature  and  time.  A,  is  the  ellipsometric  angle  for 
a  non-annealed  sample  and  the  term  <SA^“*(T„,  t„)  is  the  over  layer 
relaxation  correction.  This  correction  term  represents  the 
difference  in  A  for  the  single-film  (without  interface)  system  of 
a  non-annealed  and  annealed  overlayer  with  the  refractive  index 
calculated  from  eqns. (2-4) .  The  chosen  photon  energy  E=3.18  eV  is 
in  the  range  of  the  maximum  interface  and  the  minimum  overlayer 
sensitivity  (18) . 

Results  and  Discussion 

I.  Model  Independent  Results 


7 


Firstly,  Fig.  1  shows  that  A®*’’(E,0,T„, t,„)  is  measured  with 
enhanced  sensitivity  at  different  annealing  conditions  in  the 
spectral  energy  range  2. 8-3. 4  eV.  Figures  2  and  3  display 
experimental  results  in  terms  of  the  interface  parameter  <SAi^(T„,t„) 
defined  by  eqn. (6)  above.  Thus  Figs.  2  and  3  show  the  values  for 
SA^  measured  at  3.18  eV  versus  anneal  time  at  a  number  of 
temperatures  in  Fig.  2,  and  versus  anneal  temperature  at  two  times 
in  Fig.  3.  Other  information  relating  to  modelling  is  also 
included  in  Fig.  3,  but  discussion  of  these  items  will  be  deferred 
to  later  sections.  It  is  seen  in  Fig.  2  that  at  all  annealing 
temperatures  two  temporal  regions  of  behavior  are  present. 
Initially,  the  anneals  yield  a  fast  increase  in  SAy^  that  slows 
quickly  after  5  min  and  then  nearly  saturates.  Taking  the  data 
from  Pig.  2  at  two  anneal  times,  one  before  and  one  near 
saturation,  and  plotting  versus  anneal  temperature  reveals  a 
distinct  break  in  the  900°-970°C  range,  which  corresponds  to  the 
viscous  flow  range,  i.e.  a  temperature  above  which  viscous  flow  of 
the  oxide  is  fast  (30)  .  Thus  we  consider  that  the  viscous 
relaxation  dominates  at  the  higher  anneal  temperatures,  but  at 
lower  temperatures  we  need  to  consider  other  possible  mechanisms. 
It  is  clear  that  the  low  temperature  mechanism  accounts  for  a  large 
part  of  the  change  in  iSA;^. 

If  the  interface  region  is  treated  as  a  single  homogenous 
film,  then  the  extent  of  the  interface  is  observed  to  decrease  with 
both  annealing  time  and  temperature  as  evidenced  by  the  increasing 
Ain,.  With  both  a  decrease  in  the  interface  region  for  short  and 
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long  times  and  low  and  high  temperatures  the  models  that  are  chosen 
for  the  different  modes  of  behavior  must  show  the  same 
characteristics . 

II.  A  Working  Model  for  the  Si-SiOj  Interface. 

The  transition  from  crystalline  Si,  c-si,  to  bulk  amorphous 
SiOj  includes  at  least  two  physically  distinct  regions.  First  there 
is  a  short  range  decay  region  from  the  c-Si  in  which  the  long-range 
crystalline  order  of  c-Si  disappears,  i.e.  the  memory  of  the  Si 
order  is  lost.  The  characteristic  range  for  this  zone  is  0.1-2  nm. 
Second,  there  is  a  longer  range  decay  region  possibly  extending  to 
tens  of  nm  into  the  amorphous  SiOj  where  mechanical  stress,  density 
and  refractive  index  relax  to  the  bulk  SiOj  values.  Annealing 
decreases  the  extent  of  both  regions.  One  way  to  think  about  the 
interface  between  film  and  substrate  is  as  a  transition  region  that 
has  a  structure  with  two  major  components  each  with  a  substructure: 
the  "physical”  interface  and  the  "chemical"  interface.  The 
"physical"  interface  consists  of  a  mixture  of  substrate  and 
overlayer  constituents,  which  can  represent  microroughness  or 
protrusions  of  Si  into  the  oxide  and  even  inclusions  of  Sio, 
microcrystallites,  such  as  crystobalite  or  tridymite  (10) .  This 
layer  represents  an  order  transition.  The  "chemical"  interface 
consists  of  a  new  chemical  compound  with  a  wide  homogeneity  range, 
a  suboxide,  SiO^,  with  0<x<2.  In  the  chemical  interface  we  permit 
other  properties  such  as  stress,  density  and  refractive  index 
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(which  are  likely  related  (24)  to  relax  to  bulk  values. 

For  the  purpose  of  analyzing  the  experimental  data  in  terms  of 
physical  parameters,  we  propose  a  model  for  the  interface 
consisting  of  these  two  distinct  entities:  the  physical  and 
chemical  interfaces.  The  physical  interface  is  modeled  as  Si 
protrusions  whose  distribution  and  size  are  important  parameters. 
The  chemical  interface  is  modeled  as  a  suboxide  that  coats  the 
physical  interface.  Fig.  4  shows  a  sketch  of  our  model.  Below  we 
provide  the  parameters  that  quantify  this  working  model,  and  then 
the  experimental  data  A“’’(E,<^ ,T„, ^“*’(E,0,T„, t„)  from  Figs.  2 
and  3  are  reduced  to  yield  values  for  the  model  parameters  which 
are  compared  with  other  studies  about  the  interface  region. 
Finally,  we  discuss  the  mechanisms  that  could  diminish  both 
interfaces  with  annealing,  since  that  is  what  we  observe.  It 
should  be  understood  that  our  model  is  a  '’working"  model  based  on 
many  presently  available  research  results.  While  we  show  good 
agreement  between  the  parameters  obtained  in  this  model  and 
independent  measurements  on  the  interface  region,  we  do  not  prove 
our  model  but  only  further  demonstrate  its  reasonableness,  and  more 
importantly  provide  physically  relevant  parameter  changes  related 
to  the  condition  of  the  interface  such  as  surface  roughness  and 
stoichiometry . 

We  describe  the  crystalline  silicon  protrusions  as  hemispheres 
with  an  average  radius  R,  which  form  a  hexagonal  network  with  an 
average  distance  D  between  centers  of  the  protrusions.  Thus  the 
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protrusions  define  the  "physical"  interface  with  a  thickness  R. 
The  protrusions  and  the  region  between  them  are  covered  by  a  layer 
of  suboxide  SiO,,  0<x<2,  with  an  average  thickness  L„.  This  layer 
fonns  a  chemical  transition  zone  from  Si  to  SiOj  or  the  "chemical" 
interface.  It  is  possible  to  model  both  the  "physical"  and  the 
"chemical"  interfaces  as  one  transition  layer  with  an  effective 
thickness: 


=  (7) 

and  an  effective  dielectric  function  which  represents  a  mixture 
of  crystalline  silicon  c-Si,  silicon  suboxide  SiO,  and  the  SiOj 
overlayer  and  written  as: 

®inc  ~  ^c~Si'  ^Sio^'  ^SiO^>  ^Sio^’  ^Sio,'^  (S) 

was  calculated  using  the  BEMA,  where  the  dielectric  properties, 
c,  and  relative  volume  fractions,  f,  of  all  of  the  interfacial 
layer  constituents  are  known  a  priori.  The  dielectric  function  of 
SiO,  was  calculated  using  the  BEMA,  and  by  considering  that  SiO,  is 
a  mixture  of  amorphous  silicon,  a-Si,  and  SiO^  (38) .  The  relative 
volume  fractions  of  Si-Si  and  Si-0  bonds,  fsj^i  and  fg^^, 
respectively,  are: 


-Si-Si 


2  -X  ^  _  2x 

2  +  x  '  2+x 


(9) 


The  refractive  index  of  the  SiOj  overlayer  was  calculated  from 
eqns. (2-4) .  The  relative  volume  fractions  of  the  interfacial  layer 
constituents  f^^;,  fg^ox  fsio2  were  calculated  from  the  assumed 

geometry  of  the  interface  and  given  as: 
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where  R"  and  hj*  are,  respectively,  initial  values  for  the  average 
radius  of  the  protrusions  and  for  the  thickness  of  the  chemical 
transition  layer.  The  kinetic  coefficients  a(T„)  and  jS(T„)  are 
annealing  temperature  dependent  factors  related  to  the  rate 
limiting  processes  discussed  below.  The  powers  p  and  g  are  taken  to 
be  p=g=0.5  which  is  model  dependent  and  justified  below  based  on 
diffusion. 

The  minimization  of  the  error  function  (eqn.(5))  for  the  sets 
of  experimental  data  for  the  non-annealed  wafer  at  0=72“  and  75“ 
with  fitting  parameters  D,  R“  and  (at  x=l)  gives  the  average 
distance  between  the  centers  of  protrusions  D  =  44±4  A,  initial 
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radius  of  the  protrusions  R"  =  9.8±0.3  A  and  initial  thickness  of 
suboxide  SiO  L„  =  3.4±0.2  A.  These  results  are  in  agreement  with 
the  interface  geometry  sensitive  TEM  study  of  the  SiOj-Si,  which 
shows  that  distances  between  protrusions  at  non-annealed  interface 
are  distributed  in  the  range  40-50  A  and  heights  are  9-15  A  (1,2)  . 
Also,  photoelectron  spectroscopy  revealed  a  chemical  transition 
layer  with  a  thickness  of  2.4-4  A  (11,  44). 

A  simulation  of  the  kinetic  dependencies,  SA-^{T^,  t„)  ,  using 
egns.  (7,9-11)  with  p=g=0.5  yields  the  coefficients  ci£(T„)  and  i3(T„) 
as  functions  of  t„  or  T„  as  is  shown  in  Fig.  5.  At  annealing 
temperatures  less  then  900°C  the  kinetic  coefficient  i8(T„)  ,  which 
characterizes  the  "chemical”  interface  evolution,  is  negligibly 
small  and,  therefore,  the  "chemical"  interface  is  stable.  The 
kinetic  coefficient  a(T„)  describing  the  "physical"  interface 
evolution  or  the  silicon  protrusions  diminution,  is  saturated  and 
becomes  temperature  independent  at  T„  in  the  range  of  850-950°C.  The 
comparison  of  the  experimental  plots  of  with  simulations 
(dashed  lines  in  Fig.  2)  shows  good  agreement.  Note,  that  the 
asymptotic  form  assumed  for  the  interface  reduction  laws  given  as 
eqn. (11) ,  at  t„»l)  is  incorrect,  since  it  leads  to  a  disappearance 
of  the  interface  region  after  sufficient  time  of  annealing. 
Nevertheless,  the  phenomenological  model  satisfactorily  describes 
the  experimental  curves  of  Figs.  2  and  3  for  tens  of  minutes  and 
reflects  the  main  features  of  the  process. 

Figure  6  displays  the  decrease  of  the  extent  of  the  interface, 
i.e.  effective  thickness  from  eqn. (7)  with  a  minimum  thickness  of 
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about  3  A  realized  at  the  highest  anneal  temperatures.  The 
interface  effective  refractive  index  vs  t„  at  several  T„  is  shown 
in  Fig. 7.  Initially  both  the  index  and  thickness  of  the  interface 
layer  decrease  due  to  the  shrinking  of  the  silicon  protrusions. 
This  is  followed  by  the  slower  decrease  of  the  SiO^^  transition  layer 
which  becomes  dominant  after  considerable  reduction  of  the  height 
of  the  protrusions.  The  sharp  increase  in  the  index  seen  in  Fig. 
7  for  the  llOO^C  annealed  sample  is  indicative  of  the  index 
returning  quickly  to  the  SiO  value  as  the  fraction  of  the  Si  as 
protrusions  in  the  interface  layer  rapidly  goes  to  zero.  The  same 
effect  would  occur  for  the  other  anneal  temperatures  but  more 
slowly.  A  slower  rise  in  index  is  seen  for  the  1050"C  annealed 
sample  after  10  min  and  the  time  is  too  short  to  see  the  rise  for 
the  lower  temperature  anneals.  The  overall  shape  of  the  results  in 
Fig.  7  is  typical  of  and  dictated  by  the  BEMA  model  used  to 
interpret  the  data. 

Discussion  of  the  Interface  Model 

We  have  considered  two  interface  regimes  that  decrease  with 
annealing.  The  physical  interface  is  typified  by  the  Si 
protrusions  that  decrease  in  size  and  number  during  annealing.  The 
driving  force  for  this  is  the  minimization  of  the  interfacial  free 
energy  that  requires  an  increase  in  the  radius  of  curvature  for 
interface  protrusions  or  a  reduction  of  the  interface  area  covered 
by  the  protrusions  with  decreasing  radius,  thereby  reducing  the 
interface  width.  The  chemical  interface  is  typified  by  the 
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suboxide  and  is  reduced  by  the  conversion  of  suboxide  due  to  the 
presence  of  oxygen  which  produces  SiOj  or  the  evaporation  of 
suboxide  or  to  both,  and  by  the  relaxation  of  interfacial  stress. 
The  interfacial  stress  can  also  couple  into  the  reduction  of  the 
physical  interface  by  providing  an  additional  driving  force  for  the 
migration  of  atoms  to  reduce  the  protrusions.  We  discuss  the 
viscous  flow  and  interface  reaction  components  of  the  interface 
annealing  mechanism  separately. 

Viscous  Flow  Mechanism.  In  the  Si-SiOj  system  there  is  considerable 
evidence  that  stress  relaxation  occurs  at  high  temperature  by 
viscous  flow  (22,  24-25,  31)  with  a  viscous  flow  point  for  bulk  SiO, 
near  SSCC  (23-24,  30) .  Stress  relaxation  by  viscous  flow  decreases 
the  oxide  density  with  a  relaxation  time  r  expressed  by  eqn. (4) . 
The  characteristic  relaxation  time,  r,  is  of  the  order  of  a  few 
minutes  at  lOOO^C  and  a  few  hours  at  900°C  (22-24,  30).  This 
strongly  suggests  that  the  viscous  flow  mechanism  may  be  dominant 
at  temperatures  higher  than  -950®C,  but  becomes  negligible,  at  lower 
temperatures  and  for  annealing  times  less  then  -1  hour,  i.e.  the 
short  annealing  times  of  the  present  study.  Furthermore,  the 
temperature  activated  exponential  nature  of  a  viscous  flow 
mechanism  is  not  evident  in  Figs  3  and  5  for  temperatures  below 
950‘’C,  but  is  in  evidence  in  this  data  at  the  higher  anneal 
temperatures.  From  Fig.  3,  we  observe  that  at  T„<900°C,  (SA;^(T„)  (at 
constant  time  of  annealing)  tends  to  saturation,  and  the  kinetic 
coefficients  a(T„)  and  /3(T„)  are  saturated  at  T„=900-950°C  (Fig. 5). 
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Yet,  at  these  lower  temperatures  significant  change  also  occurs  in 
the  interface  as  evidenced  in  Figs  2  and  3 .  At  high  temperatures 
the  migration  of  Si  from  the  protrusions  occurs  via  two  driving 
forces:  free  energy  and  stress.  Another  low  temperature  mechanism 
dictates  the  interface  changes  at  low  temperatures  and  we  consider 
that  the  chemical  reactions  at  the  interface  are  likely  important. 

Interfacial  Chemical  Reaction  Mechanism.  The  reduction  of  both  the 
physical  and  chemical  interface  can  occur  by  the  reduction  in 
available  Si  at  the  interface.  Above,  the  diffusion  of  Si  atoms 
from  the  protrusions  in  response  to  the  thermodynamic  and 
mechanical  forces  was  considered.  However,  it  is  difficult  to 
understand  how  the  chemical  interface  can  be  reduced  by  this 
mechanism,  since  Si  from  SiOj  or  SiO  is  not  readily  available  to 
migrate. 

Thermal  dissociation  of  SiO^  has  received  recent  attention  (32- 
34) .  It  was  shown  that  at  elevated  temperatures  and  with  an  oxygen 
deficiency  (vacuum  or  an  inert  gas  ambient)  the  SiO^  decomposition 
takes  place  via  an  interface  reaction  (32) : 

SiOj+Si  -  2SiO 

At  temperatures  T„>-900®C  voids  are  formed  at  the  interface  after 
long-term  vacuum  annealing.  The  kinetics  of  the  void  growth  (32) 
demonstrates  that  the  oxide  decomposition  reaction  is  initiated  at 
active  defect  sites  already  present  at  the  Si/SiO^  interface. 
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Annealing  at  -750-900°C  which  is  insufficient  for  void  growth, 
transforms  these  defects  into  an  electrically  active  state  which 
seem  to  cause  low-field  dielectric  breakdown  (33,37).  The  Si 
protrusions  cannot  be  excluded  from  consideration  as  defects  that 
could  cause  the  above  SiOj  disproportionation,  since  these  sites  are 
thermodynamically  active  owing  to  the  smaller  radius  of  curvature. 
It  is  noteworthy  that  Walkup  and  Raider  (43)  did  not  observe  the 
escape  of  SiO  during  defect  formation  and  indeed  the  diffusion  of 
SiO  in  SiOj  would  likely  not  be  a  favored  process.  However,  if 
dissolved  traces  of  oxygen  or  water  are  present,  in  the  overlying 
SiOj  and/or  in  the  anneal  ambient,  this  available  oxidant  might 
remove  SiO  through  formation  of  SiOj.  In  our  experiments  we  have 
used  pure  anneal  gases  and  careful  procedures.  However,  our  anneal 
systems  are  open  flowing  systems  and  our  gases  are  from  tanks. 
Thus  trace  amounts  of  oxidants  are  present  and  can  be  dissolved  and 
accumulated  in  the  SiOj.  It  was  shown  (33)  that  with  sufficient  Oj 
(ppm  level)  present  during  inert  gas  annealing,  Oj  reoxidizes  SiO 
at  the  defects  and  prevents  the  low-field  breakdown.  This  chemistry 
may  lead  to  a  gradual  reduction  of  the  active  defect  sites  and 
causes  the  saturation  of  'SAi,^(T„,  t„)  at  T„<900°C  seen  in  Fig.  3.  From 
this  data  in  the  low  temperature  region  it  can  be  argued  that  a 
reduction  both  of  Oj  concentration  in  the  interface  layer  and  the 
active  sites  due  to  oxidation  could  be  a  rate-limiting  factor. 

In  the  anneal  temperature  region  900®<T„<970‘’C  both  mechanisms, 
viz.  relaxation  and  chemical  reaction  are  in  evidence.  We  have 
assumed  that  the  interfacial  diffusion  of  atoms  is  a  limiting 
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process  for  the  evolution  of  both  the  ”physical''  and  the  "chemical" 
interface.  At  high  temperatures  (T„>950°C)  the  "physical"  interface 
changes  by  the  diffusion  of  Si  atoms  away  from  protrusions  and 
viscous  flow  of  SiO,  (or  SiOj)  in  the  opposite  direction.  At  1000°C 
the  viscous  flow  relaxation  time  is  a  few  minutes,  but  as  seen  in 
Fig. 2  the  interface  continues  to  change  for  longer  times  with  a 
slower  rate.  After  the  fast  and  ample  reduction  of  the  extent  of 
the  physical  interface  for  t„<~10  min,  the  interface  evolves  slowly 
(with  a  «  a)  by  the  reduction  of  the  chemical  interface  via  the 
interfacial  diffusion-limited  reactions  of  decomposition  and 
oxidation  of  suboxide.  At  moderate  temperatures  (T„<900°C)  the 
decrease  of  SiO,  (x<2)  is  diffusion  limited  (36)  likely  by  the 
diffusion  of  oxidant  through  the  silica  network  to  the  interface. 
Thus,  it  was  reasonable  to  use  a  diffusional  square-root  time  law 
for  the  "physical"  and  "chemical"  interface  evolution,  i.e. 
p=g=0.5,  and  a  good  fit  was  found.  According  to  the  model 
developed  here,  annealing  in  an  oxidant  free  atmosphere  at  low 
temperatures,  where  the  reoxidation  reaction  2Si0+02  -♦  2Si02  cannot 
take  place  will  result  in  no  annihilation  of  the  electrically 
active  sites. 

Conclusions 

A  novel  enhanced  interface  sensitivity  immersion  spectroscopic 
ellipsometry  technique  was  applied  to  study  the  Si/Si02  interface 
annealing.  Two  distinct  stages  of  interface  evolution  were  found: 
a  fast  initial  stage  with  interface  microroughness  reduction 
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followed  by  a  slow  decay  of  the  interfacial  suboxide.  The 
interfacial  suboxide  reoxidation  reaction  dominates  the  interface 
evolution  in  the  moderate  temperatures  region  (T„  <  900°C)  and 
viscous  flow  becomes  dominant  at  elevated  temperatures  (T„  >  950°C)  . 
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FIGURE  CAPTIONS 


Figure  1.  The  experimental  dependence  of  the  ellipsometric  angle, 
A,  on  photon  energy,  E,  at  different  annealing  times. 

Figure  2.  The  experimental  dependence  of  the  effective  interface 
parameter  ,  A-^,  on  annealing  time  at  a  number  of  annealing 

temperatures.  A  simulated  dependency  is  shown  with  the  dashed 
curve . 

Figure  3.  The  experimental  dependence  of  the  effective  interface 
parameter,  A^,  on  annealing  temperature  for  two  annealing  times. 

Figure  4.  Schematic  representation  of  our  interface  model. 

Figure  5.  The  dependence  of  the  kinetic  coefficients  a  and  S  on 
annealing  temperature. 

Figure  6.  The  dependence  of  the  effective  interface  thickness  on 
annealing  time  at  a  number  of  annealing  temperatures. 

Figure  7.  The  dependency  of  the  effective  interface  refractive 
index  on  annealing  time  at  a  number  of  annealing  temperatures. 
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Protrusions 


